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We give an update on our calculation of Bk using HYP-smeared valence staggered quarks. We 
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1. Introduction 

An accurate result for the kaon B-parameter is important both for its phenomenological im- 
pact ^ and as a bellwether of the success in incorporating chiral symmetry and controlling sys- 
tematics in lattice calculations. Calculations using several different fermion methods are underway, 
with the present best result using domain- wall fermions We are pursuing a calculation using 
improved staggered fermions. This has the advantage of being computationally cheap, but the chal- 
lenge of deahng with the effects of taste-breaking in a context where there is operator mixing 

We use the standard staggered action with HYP-smeared links for our valence fermions, which 
reduces taste-breaking by a factor of 3 compared to asqtad quarks |Q ^. We use the MILC lattices 
generated with 2 + 1 flavors of asqtad sea quarks. For the coarse MILC lattices, on which we focus 
here, the resultant taste-breaking in the pion masses, while reduced compared to asqtad quarks, 
remains large enough that we must use the standard power-counting of staggered chiral perturbation 
theory, in which p^. The complications that this introduces have been explained in Ref. 

Our method for calculating Bk using wall sources is explained in Ref. For each lattice 
we use 10 valence quark masses running from Ki mf^^^ down to « mf^^^ / 10 in equal steps, and we 
calculate Bk and tuk for the Goldstone taste for all 55 quark-mass combinations. Results for one- 
loop matching factors using the mixed action are not yet complete so we use tree-level matching. 
For this, and other reasons to be explained, all results obtained here should be regarded as very 
preliminary. They are essentially our first pass at fitting the data, which we are using to inform 
subsequent fitting and to determine where improvements in statistics are needed. 

2. Staggered chiral perturbation theory for a mixed action 

Staggered chiral perturbation theory (S;^PT) ||7|, ||] incorporates discretization errors into the 
chiral expansion, and in particular includes the effects of taste-breaking and rooting. For our light- 
est kaons, such effects are comparable to those coming from the explicit chiral symmetry breaking 
due to quark masses, and thus enter S;tPT at LO. A major effect in S;tPT is that chiral loops 
(which begin at next-to-leading-order [NLO]) must be evaluated with the masses of the pions of 
the appropriate tastes, rather than a common mass. For Bk, which involves an insertion of the weak 
Hamiltonian, one must also deal with mixing between operators having different tastes. This leads 
to a significantly larger number of unknown coefficients multiplying NLO terms than are present 
in the continuum, as explained in Ref. [^. 

The S;^PT analysis of Ref. [||] does not apply directly to our set-up, however, because we 
use a mixed action. Following the methods developed in other mixed-action contexts |^, we have 
generalized the results of Ref. to our setup. This turns out to be straightforward. Here we only 
give a summary — further details will be presented in Ref. [[To|]. 

There are three classes of effects resulting from using a mixed action. The first, which cor- 
responds to the short-distance parts of sea-quark loops, is simply that the coefficients multiplying 
various terms in the S;t;PT expression for Bk will change. Since these coefficients were previously 
unknown, however, this has no practical impact — one trades one set of unknown coefficients for 
another. 
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The second class comes from loop diagrams involving mixed pions — those composed, say, of 
a valence (HYP) quark and a sea (asqtad) antiquark. These, however, are absent for Bk at NLO. 

The third class involves loop diagrams with pions composed of sea-quarks alone. For Bk at 
NLO these loops all include "hairpin" vertices, for this is the only way in which the long-distance 
part of sea-quark loops can enter. The effect of using a mixed action then boils down to the need to 
distinguish between three types of hairpin vertices — valence- valence, valence-sea and sea-sea — all 
of which would be the same if the same valence and sea quarks were being used. It turns out that 
this has an impact only for tastes V and A (the taste singlet hairpin coupling to a particle which 
is being integrated out anyway), so that one ends up with 6 hairpin parameters (3 for each of two 
tastes) instead of 2. We name the hairpin parameters 5^'^', and 5^", with B = Aor V, and the 
superscript indicating the types of quark involved. They appear at NLO in the combinations 

= iS's -'f/S's'' , = - {5'^'] /5s'' • (2.1) 

These only enter in the expressions for non-degenerate quarks — those for degenerate quarks in the 
kaon do not involve hairpin vertices. Specifically, in the contribution denoted ^^^^ in Ref. [^], and 
given in eq. (50) of that paper, one must, for both B = V and A, make the substitution 

5^-5r\ (2.2) 

and add the following new term 

a^Sr'^^j^ (^2^M^ + +£(F,)) . (2.3) 

Here we use the notation Xb (Yb) for the mass-squared of the flavor non-singlet valence pion with 
taste B and composition xx (yy), where the kaon itself has the composition xy (and taste P). In the 
notation of Ref. Xb = and Yb = my^. The functions i and I are chiral logarithms, and are 
defined in Ref. 

The original hairpins, enter through their (unchanged) contributions to the masses of the 
flavor-singlet mesons, r\B and T]g. 

The hairpin parameters are a measure of taste-symmetry breaking, and we expect them to sat- 
isfy a similar hierarchy to that we observe in the pion spectrum, namely 8b^' / 8b' ~ 1/3. In words, 
hairpins correspond to quark-antiquark pairs communicating through intermediate gluons, and the 
taste-breaking component is, by construction, reduced for our HYP-smeared valence quarks. Based 
on this argument we also expect that (5g^'V ~ 1/3. Combining these expectations we find 

and S^^^^O. (2.4) 

If these expectations are accurate, then using a mixed action has essentially no impact on chiral 
fitting, because the new terms proportional to 8b [eq- (2.3)] can be dropped, and the terms 



proportional to 8b were present anyway with an unknown coefficient. 
3. Fitting strategy 

The NLO expression for Bk in partially-quenched mixed-action S;^PT takes the form 

16 

BK = Y.Cifi, (3.1) 

(=1 
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in which ft are known functions and c, are coefficients to be determined. Of the 16 unknown 
coefficients, 1 (cq) appears at LO (and is the value of Bk when a —^0 and then the chiral limit is 
taken), 4 are the NLO low-energy constants (LECs) present in the continuum, and the remaining 1 1 
are LECs due to lattice artifacts and truncated perturbative matching. For degenerate quarks these 
numbers reduce to 9 coefficients = 1 LO + 3 continuum LECs + 5 "lattice LECs". This counting 
is for a single lattice spacing — the dependence of the artifacts on a involves a mix of a^, and 

The functions /)■ depend on the masses of the flavor-non-singlet valence pions of all tastes and 
all compositions (xx, yy and xy), which we determine as part of our calculation [^. In addition, we 
need the masses of the sea-quark i£ and ss flavor non-singlet pions with taste V, A and /, and of the 
flavor singlet rj^ and Tj^ for 6 = V and A. These we take from the results of the MILC collaboration 
(including the axial and vector hairpin vertices 5^^) [11]. The final inputs we need are a (we use 



the MILC values) and / (which, for the moment, we simply set to / = 132MeV). 

Although we have 10 degenerate and 45 non-degenerate data points on each lattice, a direct fit 



to eq. ( p.ll ) is difficult, since many of the fit functions are similar. As a first stage, therefore, we use 
only a few representative "lattice" contributions, while keeping all the continuum terms. We also 
use Bayesian priors to constrain some terms, since we know their orders of magnitude. 
The fit functions we use are, firstly, 

/l = 1 + ^ i^^^onn + J^discr''''''''"'^ , (3.2) 

where G = m^.p, ^^cotm,disc come from quark connected/disconnected 1-loop diagrams and are 
given in Ref. [^, and the superscript indicates that we keep only the contributions from chiral 
operators present in the continuum. We do, however, include taste breaking in the pion masses in 
these contributions. We set the scale in the chiral logarithm to /i = 1 GeV. Next we include the 
continuum analytic terms (with the ;^PT scale set to A = 1 GeV): 

fi = G/A^ , /3 = {G/AY , fs = {Xp - Ypf/{GA^) , 

Note that /s is a NNLO contribution, but is needed to fit our data up to the highest quark masses. 

Finally, we include three representative discretization terms. We use the contribution to ^conn 
containing taste-J pions, which contains two parts with independent coefficients [^: 

/6 = ^ {l{XT)^i{YT)-U(m\j)) (3.3) 

and the remainder which we call /4. We also include the likely dominant contribution to ^disc, 
which is that proportional to 5^^\ We set c? = {2Cf +CY)a^5^'^^ /{n^f), and /? is the co- 
efficient of this term in eq. (50) of Ref. Note that /s^e.v contribute only for non-degenerate 
quarks. 

Based on the power-counting of Ref. [^, the coefficients should have the magnitudes: 

Ci «C2 f«C3 «C5 ~ ^(1) 

C4 R:iC6 ?aAgcz)(«^ecD)^ = 0.003 to AgcD«s = 0-01 GcV^ on the coarsc lattices. (3.4) 
C7 AgcD(«^ecz))^ = 0.0003 to AgcD«.s^ = 0-001 GeV^ on the coarse lattices. 
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Figure 1: Fits of (with degenerate quarks) vs. mj^ on the MILC coarse lattices with ani( = 0.01 and 
flffjs = 0.05. The fits are to continuum PQ;)fPT (left) and SjPT (right), as described in the text. 

Here, we use Aqcd ~ 0.3 GeV and = cc^{lJ. = l.6GeV) w 0.36. The smallness of C4 6,7 is 
somewhat offset by the fact that /4,6,7 are logarithmically divergent in the chiral limit. 

4. Examples of fits 

We first try fitting to the degenerate data without the lattice terms, so that only ci_3 are non- 
zero. We compare using only the Goldstone-kaon mass in loops ("PQ;^PT" fit) to using the appro- 
priate combination of all tastes ("S;tPT"). The resulting fits are shown in Fig. The S;^PT fit is 
significantly better, because its smaller curvature at small quark masses more accurately represents 
our data. This is the familiar "softening" of the chiral logarithms caused by the heavier masses of 
non-Goldstone taste pions. 



fit type 


c\ 


C2 


C3 


C4 


C5 


C6 




X^/dof 


D-T4 


0.39(1) 


-.01(1) 


0.89(17) 










.05(11) 


D-BT4 


0.390(3) 


-.011(7) 


0.90(4) 


.0001(4) 








.07(16) 


ND-BT4 


0.390(2) 


-.011(1) 


0.90(2) 


.0001(2) 


.13(7) 


-.011(7) 


.0018(8) 


.06(7) 


ND-T2 


0.31(7) 


.6(5) 


0.23(50) 


.002(2) 


.04(4) 


-.005(4) 


.0010(5) 


.03(2) 



Table 1: Fitting parameters. Fits are described in the text. 



The degenerate data itself shows no indication of a logarithmically divergent contribution (as 
would be produced by /4). In a first attempt to quantify this, we try and make use of the fact that 
the dominant contribution from f^. is for the lightest few kaon masses. Thus we drop the lightest 
mass point from the S;^PT degenerate fit, giving a fit we name "D-T4", whose parameters we list 
in Table [l[ They have the expected magnitudes (except perhaps for C2, but this is scale dependent, 

'We fit using an uncorrelated with errors in fit parameters determined by jackknife. Thus only the relative 
goodness of fit can be estimated, but not the absolute goodness of fit. 
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Figure 2: ABk vs. m|-, including degenerate (red) and non-degenerate (blue) data. 

and becomes 0.4 if /i = 0.77 GeV instead of /i = 1 GeV). We then do a fit to all 10 degenerate 
points including /4, but with ci_3 constrained, using the values and errors from D-T4 as Bayesian 
priors [|l|]. This fit ("D-BT4" in Table [l|), has a very small C4, with the /4 term making no more 
than a 1% contribution. Finally, we fit to the full data set (55 points) now including fs-j, but with 
ci_4 constrained using the results from fit D-BT4. The resulting fit we call "ND-BT4". We have 
also done an unconstrained fit using all 7 to all 55 data points — fit "ND-T2" in the Table. 

The difference between fits ND-BT4 and ND-T2 indicates the size of the present uncertainty in 
the coefficients. While this is substantial, it is encouraging that the coefficients have sizes roughly 



consistent with the estimates (3.4). Furthermore, when we use the fit form to determine our best 



estimate for the continuum Bk for physical kaon masses, we obtain consistent values: 0.76(6) and 
0.67(4).-' It is also encouraging that the overall fits looks reasonable. We illustrate this in Fig. 
where we compare the residuals, ABk{x) = Bk{x) — f{x) for a continuum PQ;^PT fit (c4,6,7 = 0, 
Goldstone kaon masses only) to the S;^PT fit (ND-BT4). We notice a significant improvement in 
the fitting quality using S;tPT. 

We have repeated this analysis on four other coarse MILC ensembles and one fine ensemble. 
Some results are collected in Table ^ All we can conclude at this stage is that there is a rough 
consistency between different coarse ensembles, and that discretization errors are not enormous. 
We stress again that these results use tree-level matching and so are very preliminary. 

5. Conclusion 



We have taken the first stab at fitting our mixed-action results for Bk- We clearly have a 
lot of work to do to control the systematic errors (and at this stage cannot quote a continuum 
result with all errors estimated). It is important to keep in mind, however, that the main goal of 



^The relatively large size of c-j in fit ND-BT4 is possible, in part, because of a cancellation with the eg contribution, 
and it may be better to directly constrain both these coefficients to have smaller magnitudes. 

■'These values are obtained by setting C4 g 7 = 0. We stress, however, that they still contain those discretization 
errors that are absorbed into cj , since they are based on fits at a single lattice spacing. 
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a (fm) 


anil /am., 


geometry 


ens 


Bk( tree, ND-BT4) 


Bui tree, ND-T2) 


0.12 


0.03/0.05 


20^ X 64 


564 


0.83(6) 


0.63(4) 


0.12 


0.02/0.05 


20^ X 64 


486 


0.72(6) 


0.71(4) 


0.12 


0.01/0.05 


20^ X 64 


671 


0.76(7) 


0.67(4) 


0.12 


0.007/0.05 


20^ X 64 


651 


0.87(5) 


0.62(3) 


0.12 


0.005/0.05 


24^ x64 


509 


0.80(4) 


0.66(3) 


0.09 


0.0062/0.031 


28^ X 96 


995 


0.72(4) 


0.62(3) 



Table 2: Comparison of tree-level from different ensembles (preliminary) 



the fitting is to provide a reasonable extrapolation formula to the physical quark masses. We are 
less interested in the coefficients themselves (except for cq, which can be compared to results from 
large Nc approaches). We are also investigating fits based on (mixed-action staggered) SU (2) chiral 
perturbation theory, which appear to be much simplified because discretization terms are of NNLO. 
We also expect that the more extensive data on the fine lattices which we are presently collecting 
should be more straightforward to fit. Finally, we are improving our statistics on several of the 
coarse and fine ensembles. 
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